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Multifractal scaling of critical wave functions at a disorder- driven (Anderson) localization transi- 
tion is modified near boundaries of a sample. Here this effect is studied for the example of the spin 
quantum Hall plateau transition using the supersymmetry technique for disorder averaging. Upon 
mapping of the spin quantum Hall transition to the classical percolation problem with reflecting 
boundaries, a number of multifractal exponents governing wave function scaling near a boundary 
are obtained exactly. Moreover, additional exact boundary scaling exponents of the localization 
problem are extracted, and the problem is analyzed in other geometries. 



I. INTRODUCTION 



In a system of non-interacting electrons disorder can 
induce transitions between -metallic (delocalized) and 
insulating focalized) phases! The scaling theory of 
localization^ predicts that in two dimensions (2D), all 
single particle electronic states are typically localized by 
arbitrarily weak disorder and consequently, there are no 
metal-insulator transitions for non-interacting electrons 
(except in the situations discussed below). This scal- 
ing idea can be substantiated quantitatively through a 
non-linear sigma model (NLSM) formulation of_the prob- 
lem using either the replicao or supersymmetryQ method. 
The NLSM approach has been very successful in provid- 
ing a quantitative understanding of the metal-insulator 
critical point in dimensions greater than two (for a re- 
view, see H 

There are two well known exceptions to the above dis- 
cussion. One is the 2D integer quantum Hall (IQH) 
plateau transition and the other is the metal-insulator 
transition in 2D systems with spin-orbit scattering (the 
so-called 'symplectic symmetry class'). These possess 
phases of both localized as well as extended states in 
2D. On one hand, both these transitions can be for- 
mulated as NLSM field theoriesftEl (containing an ad- 
ditional, topological term in the IQH case). But, on 
the other hand, these field theories cannot be used to 
perform analytical (or even perturbative) calculations of 
critical exponents or other (universal) physical quanti- 
ties at the transition. This is because the long-distance 
physics relevant for the properties at the transition is 
governed by the strong-coupling regime of these theories, 
which is in general not tractable. Indeed, such a cal- 
culation would require a non-perturbativc approach to 
these problems. Some progress on these types of prob- 
lems has been possible in recent years due to the appear- 
ance of models which are analytically tractable even at 
strong coupling and which exhibit similar localization- 
delocalization (LD) transitions, but belong to symme- 
try (and universality) classes different from the ones dis- 



cussed above. Indeed, semijial work due to Zirnbauer, 
and Altland and ZirnbaueitJiI!J, which appeared a little 
more than a decade ago, has identified a total of ten 
symmetry classes, describing in principle all possible be- 
haviors of non-interacting quantum mechanical particles 
subject to random disorder potentials. Examples of such 
analytically tractable 2D localization-delocalization (LD) 
transitions iiplude 2D Dirac fermions subject to. 
dom AbcliantiJ and non-Abelian gauge potentj^lsE 



the spin quantum Hall (SQH) plateau transitionliiHi^El 
(reviewed briefly in Section [ij below), and-problems with 
a special so-called 'sublatticc symmctryliS. 

LD transitions are different in several respects from 
continuous phase transitions occurring in non-random 
systems, such as the ordering transition in a clean fer- 
romagnet. One difference is the absence of an obvious 
order parameter in many cases. The disorder-averaged 
density of states which would be expected to play the role 
of an order parameter, is not criticafl at many conven- 
tional LD transitions, including those mentioned at the 
beginning of this Section. Another very important differ- 
ence is that the wave functions at the transition .e xhibi t 
what is known as multifractal scaling bchavionj'tj'ErjEil. 
This means that the disorder average of the q-th. power 
of the (square of) the critical wave function amplitude 
scales with the system size L with a non-trivial expo- 
nent that has a non-linear dependence on the power q. 
In field-theoretic language, multifractality can be under- 
stood as the existence of a certain infinite set of fields 
in the theory which exhibit infinitely many independent 
anomalous dimensions with certain very special ('convex- 
ity') propcrtiesc2l. In contrast, analogous infinite sets of 
fields cannot exists in 'conventional', non-random phase 
transitions. (In 2D, for example, such clean phase tran- 
sitions are known to possess onlyj-a finite number of in- 
dependent anomalous dimensionsE3.) On one hand, the 
multifractal scaling properties of wave functions at the 
2D IQH plateau transition and at the 2D symplectic class 
metal-insulator transiting iiave been calculated numeri- 
cally to high precisioi£J'E3'E3. On the other hand, the en- 
tire multifractal spectrum has been obtained analytically 
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for Dirac fermions in random gauge potcntialaiiltJO 
while the first few multifractal wave function moments 
have been calculated analytically for the spin quantum 
Hall transitionBEl. 

The aforementioned properties, while being important 
characteristics of the LD critical points, do not exhaust 
all their important universal physics. It is well known 
that even for non-random, conventional phase transi- 
tions, the geometry of the system has a crucial effect 
on the critical behaviorEil. Critical exponents as well 
as correlation functions depend on the overall geometry 
and the location where these quantities are measured. 
In particular, correlation functions of quantities located 
near the boundary of the system exhibit in general dif- 
ferent scaling behavioifijOO as compared to bulk cor- 
relation functions (and this depends in general also on 
the imposed boundary conditions). These aspects can 
be studied in great detail especially in the casfi-pf criti- 
cal 2D systems which are conformally invariant^. These 
ideas thus emphasize the importance of boundary effects 
in LD transitions, particularly in 2D where these tran- 
sitions are believed to be in general described by con- 
formally invariant fixed points (see also Ref. ^6|). The 
SQH transition, being analytically tractable, provides us 
with a playground where we can understand the effect 
of both multifractality and boundary criticality at LD 
transitions. In the process we are led to a more general 
framework where one considers multifractal behavior of 
wave functions near boundaries, giving rise to the notion 
of 'surface (boundary) multifractality'. The power of con- 
formal invariance enables us to extend this idea to more 
complicated geometrical settings. In the present paper 
we derive several exact exponents characterizing these 
properties for the SQH transition. 

To summarize our main results, we note that 
within a second-quantized formulation of the conven- 
tional Chalker-Coddington network model for the SQH 
transitionEjO, a reflecting boundary (where the system 
simply ends) preserves the full sl(2|l) supersymmmetry 
present in the, disorder-averaged bulk problem. After 
the mappingiH to bond percolation this amounts to a 
reflecting boundary condition for the perimeters (hulls) 
of the percolation clusters. Here wc express correlation 
functions characterizing wave function multifractality in 
terms of the generators of the underlying sl(2|l) super- 
symmmetry algebra. Using the representation theory of 
this superalgebra together with the percolation mapping 
we are able to convert the correlation functions to clas- 
sical percolation probabilities. Using previously known 
scaling exponents for percolation we are able to calcu- 
late two non-trivial multifractal exponents, convention- 
ally denoted by A2 and A3, both on the boundary and in 
the bulk, within the supersymmetry approach developed 
in Ref. [I?]. (The values for A2 and A3, in the bulk had 
been obtained in Refs. ^ and [30] using methods very 
different from the supersymmetry methods used in the 
present paper.) The inability to calculate higher multi- 
fractal moments is explained within our (supersymmetry) 



formulation. The scaling behavior of the averaged local 
density of states and point-contact conductance are also 
found at the boundary. Using suitable conformal invari- 
ance arguments, we are able to translate the boundary 
('surface') exponents to a wedge geometry and find the 
corresponding corner exponents (see Section VI for the 
various exact values). 

The organization of this paper is as follows. In Section 
H we briefly present the physics of the SQH effect. In Sec- 
tion H|, we discuss the supersymmetry (SUSY) formula- 
tion of the Chalker-Coddington network model for the 
SQH transition in the presence of boundaries. In Section 



IV, we elucidate how supersymmetry in the SQH prob- 
lem can be used to study (low-order) multifractal wave 
function moments, both near boundaries and in the bulk. 
(As mentioned above, exponents of low-order multifrac- 
tal moments in the bulk have been obtained previously 
in Ref. ^ using very different methods.) We also ex- 
plain, using the supersymmetry approach, following Ref. 
p"7| , the technical reasons why this computation can be 
done only for low multifractal moments (up to the third 
power of the wave function intensity) and not for the 
whole multifractal spectrum. The reasons appear to be 
different from, and complementary to the ones discussed 
111 Ref. |0[ In Section Q we discuss the scaling behavior 
of the local density of states and of the point contact con- 
ductance near a boundary. In Section VI , we extend our 



discussion to more complicated geometries. In the final 
Section, we conclude by discussing certain implications 
of our results and a number of remaining open issues. 
Some of the results derived in the present paper were 
announced in the Letter Ref. B6l 



II. PHYSICS OF THE SPIN QUANTUM HALL 
TRANSITION 



The spin quantum Hall (SQ H) transition was first 



studied, numerically, in Ref. 15, A simple physical pic- 
ture of the SQH effect was given in Ref. [ll| In this Sec- 
tion we briefly review, for completeness, the basic physics 
of the SQH effect. 

Let us consider a lattice version of the BCS Hamilto- 
nian describing a singlet superconductor (classes C or CI 
in the notation of Refs. |9] and fio| ): 



n 



ij 



The first term here describes (in a second-quantized lan- 
guage) hopping of electrons of either spin a =1,1 be- 
tween lattice sites i and j, as well as possible on-site 
potentials (for i = j). The other two terms describe BCS 
singlet pairing. Provided that the gap function satisfies 
Ajj = Aji, the Hamiltonian (|l|) is SU(2) invariant and 
commutes with the three generators of the global SU(2) 
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spin rotations (total spin): 

n 



(2) 



Being the Hamiltonian of a superconductor, Eq. ([[]) 
does not conserve the particle number, or charge. Phys- 
ically this is due to the existence of the pair condensate 
which may exchange electrons into holes in processes sim- 
ilar to Andreev reflection. Not being a conserved quan- 
tum number, the quasiparticle charge cannot be trans- 
ported by diffusion. Consequently, there is no notion of 
electrical (charge-) conductivity. However, the spin of 
each quasiparticle is conserved by the Hamiltonian (Q), 
and one may define the spin conductivity as the linear re- 
sponse coefficient between the (say) z-component of the 
spin current and the gradient of a Zeeman field along the 
z direction: 



Ji 



(3) 



Here g is the gyromagnctic ratio, and \ib is the Bohr 
magneton. 

When the gap function A is complex, time-reversal 
symmetry is broken, and then the Hall (transverse) spin 
conductivity a xy may be non- vanishing. It is convenient 
to perform a particle-hole transformation on the down- 
spin particles: 



(4) 



Under this transformation the z component of the gen- 
erator of spin rotations (||) becomes simply the total 
number of d-particles, | Yli a ^la^ia- Thus the trans- 
formation interchanges the role of particle number and z 
component of spin, with the result that the transformed 
Hamiltonian conserves the number of d particles: 



'3/3- 



(5) 



Here Hij is the so called Bogoliubov-de Gennes (BdG) 
Hamiltonian, a 2 x 2 matrix in the spin space: 



/ ty Ajj 

iJ «i — \ a* _+* 



(6) 



Since the e?-particlc number is conserved we can use 
a single-particle description. The external non-random 
Zeeman field B z in the z direction, in this alternative de- 
scription maps onto a simple shift e in the Fermi energy 
as 



9VbB z 



(7) 



The spin conductivity in the original problem becomes 
the usual electrical (charge) conductivity of the d parti- 
cles. Note that the role of the elementary charge is now 
played by ft/2, so that the natural quantum of the spin 



conductivity is (h/2) 2 /2ttH = h/8n. In what follows we 
use units in which h = 1. 

Invariancc of the Hamiltonian in Eq. (|l]) under the 
global SU(2) symmetry implies the following property of 



the matrix Hij in Eq. (^]), 



(TyHijC 



-HI 



(8) 



which has important consequences. Firstly, it implies 
that the eigenvalues of H^ always appear in pairs as 
(e, — e). In other words, the single particle energy spec- 
trum is symmetric about e = 0. The ground state of 
the problem is obtained by filling up all the negative en- 
ergy (d-) single particle states. The positive energy (d-) 
particle and hole excitations on top of this ground state 
are doublets under the global SU(2) symmetry. A Zee- 
man term, or non-zero energy shift (see Eq. (Q)), breaks 
the global SU(2) symmetry and splits the degeneracy be- 
tween the above-mentioned (d-) particle and hole states. 
Secondly, Eq. (^) implies a relation between the retarded 
and advanced Green's functions of the BdG Hamiltonian 
or the corresponding network model, (see Eq. ( |l6|) be- 
low). 

A clean (that is non-random) d x i_ y i + id xy -wave 
superconductor (which breaks time-reversal symmetry) 
can, for example, be obtained on the lattice with a gap 
function Ay whose Fourier transform is (square lattice 
"d-wave"): 



Aq (cos k x — cos k y ) 



iA xy sin k x sin k y , 



(9) 



It is easy to seeB that this represents at low energies^ 
(2+l)-dimensional Dirac fermion of mass A xy , leadingO 
to a quantized spin Hall conductivity (in units of h/8ir): 



< y = 2sgnA a 



(10) 



This is consistent with the existence of two spin-current 
carrying states at the edge of a system with a boundary 
(say, to vacuum). 

These edge states, being chiral, survive the addition 
of weak disorder, so the quantization of a xy persists in a 
disordered d-wave superconductor with broken time re- 
versal invariance, at least for weak disorder. In general, 
there can be two possible localized phases — spin in- 
sulator, and SQH phase, — distinguished topologically 
by the quantized value of the spin Hall conductivity. A 
transition between them, the generic SQH transition, is a 
localization-dclocalization transition similar to the IQH 
plateau transition, where a xy jumpa by two units. It 
is a quantum percolation transitiontll of the edge states 
formingEJ at the interfaces between two topologically dis- 
tinct regions. 

There is one essential difference between the SQH and 
the IQH case. In the latter case, the mean single-particle 
density of states is non-vanishing on either side and at 
the transition. In the SQH case, on the other hand, the 
density of states of e?-particles, p(e), .-vanishes at zero en- 
ergy, e = 0, in both localized phasesL3 as e 2 ; right at the 
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transition, that is at e = 0, it vanishes (as might have 
been expected) with a non-trivial exponent: 



(11) 



The exponent a was calculated exactly in Ref. [l?] to be 
1/7. In Section |y| below we will present a derivation of 
this result, as well as its analogue for the local density of 
states at a boundary of the SQH system. 



III. SQH NETWORK MODEL 
A. Description of the network 

Network models have been very convenient for both 
numerical and analytical work on various disordered non- 
interacting fermion problems. The best known such, 
model is the Chalker-Coddington (CC) network modelt£l 
for the integer quantum Hall plateau transition. 

Similar network models can be constructed for other 
localization problpns- including the chiral metalocd, the 
SQH transitionHtSllZl, the random-hmnd Tsing model and 
the thermal quantum Hall effectcj'oao. 

The SQH network consists of a lattice of directed links 
and two types of nodes, A and B, forming a square lattice 
(see Fig. [l) on which spin-1/2 particles at energy e = 
can propagate. Uni-directional propagation through each 
link is represented by a random SU(2) matrix. As in the 
case of CC network, to study the critical behavior at 
the SQH transition it is sufficient to introduce disorder 
only for propagation along the links, while all the nodes 
can be taken to have the same (non-random) scattering 
matrices. The node scattering matrices are diagonal in 
the spin indices: Ss = Ssi <8> Ssi, 



(1 



- t 2 
— tSo 



1/2 



tSa 

(i - 4, 



,1/2 



(12) 



where S = A,B labels whether the node is on the A or 
the B sublattice, and a =\, \ labels the spin-index of the 
propagating particle. Apart from the case of boundary 
nodes (which will be treated in detail later), the remain- 
ing network is isotropic (invariant under 90 degree ro- 
tation of the lattice) when the scattering amplitudes on 
the two sublattice nodes are related by t 2 Aa + t 2 B(J = 1 . 
The critical point of the isotropic network is located 
at tAo = tsa = l/v2- Varying tga while keeping 
t A(J + t Ba = 1 and ts] = tgi drives the system between 
a spin insulator and a SQH state and the spin Hall con- 
ductance jumps from to 2. Taking ts^ ^ t$i breaks 
the global SU(2) symmetry and splits the transition into 
two ordinary IQH transitions each in the unitary symme- 
try class (investigated, for example, in Ref. |39| ). In this 
paper we consider only the spin-rotation invariant case 
with t S i = tsi = t s - 




FIG. 1: SQH network with a vertical reflecting boundary. 
The black squares represent the non-random node scattering 
matrices [Eq. (113)]. The black circles on the links represent 
the random SU(2) scattering matrices. 



B. Green's functions and symmetries 

Network models can be studied using either a first 
quantized or a second quantized formalism. The first 
quantized method is adopted in Refs. [l6| and It is 
also very useful for numerical work. Here we use it to 
derive certain symmetry relations. The rest of the pa- 
per will employ the second quantized formalism (see next 
Section). 

All physical quantities of interest in our problem such 
as wave function correlators and conductance can be ex- 
pressed in terms of the Green's functions. Usually, the 
latter are represented in continuous time notation, but 
the network models use the discrete time analogs. In par- 
ticular, network models use a single step time-evolution 
operator hi which acts on the single particle wave func- 
tion ^(r, t) at discrete time t to give the wave function 
at time t + 1, where r denotes a link of the network. For 
a network with N links, hi is (due to the spin- index) a 
2N x 2N matrix, which represents a finite-time version 
of the infinitesimal time-evolution operator. Thus U = 
cxp{i7i} where the Hamiltonian TL (which describes the 
time evolution of the edge states separating 'puddles' of 
topologically distinct regions) has the same symmetry 
properties as the underlying BdG Hamiltonian in Eq. 
(g). Therefore, due to the property in Eq. (||), U is 
a (unitary) symplcctic matrix and as such satisfies the 
condition 

U- 1 = a y U T a y . (13) 

where a y is the conventional Pauli matrix. 

Now we can write the retarded (advanced) Green's 
function Gn(ri,r-2) {Ga{t\, ^2)) as the resolvent of the 
operator hi : 



[Gfl(ri,r2)] a/3 
[G A {rur 2 )\ ap 



1 



1 - zU 

1 



1 - z^U 



T2,(3 



(14) 
(15) 
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Here z = e^ €+ * 7 ^ where e is the energy and 7 is a finite 
level broadening. Since the SQH transition of interest to 
us occurs only at zero energy, we will set e = and hence 
z = e 1 from here on. In the above equation, r\,r2 
and a, (3 denote the network link and the spin index, 
respectively. 

Making use of Eq. (pi) in Eqs. (Ill ||) leads to 



the following relationship between advanced and retarded 
Green's functions 

[G A (r 1 ,r 2 )] a ,- 5 a > a S rur2 = e Q/3 [G R (r 2 , n)]^^^, 

(16) 

where e a p is the antisymmetric Levi-Civita tensor. The 
above equation will turn out to be crucial for our calcula- 
tions since it implies that we need not introduce separate 
operators for the advanced sector Green's functions but 
can work solely with retarded ones (as we will sec below) . 



C. Second-quantized description in the bulk 

All network models can also Jae— studied using the 
second-quantized SUSY techniqudid'cZI. 

The second-quantization of the SQH network is de- 
cribed in Rcf. [Tt] and wc review the main steps here. The 
basic idea is quite similar to the transfer matrix formula- 
tion of various 2D statistical mechanics problems such as 
the Ising model (for a review, see Ref. ^). The presence 
of disorder introduces a number of additional features 
which arc sketched below. 

All physical quantities in the network model includ- 
ing conductance, wave function amplitudes etc. may 
be expressed in first quantization in terms of sums over 
paths on the network. Such a sum may be written in 
second-quantized SUSY language as a correlation func- 
tion, (...) = STr[ . . . {U bU a) Lt \ where the supertrace 
'STY' contains the row-to-row transfer matrices Ua and 
Ub and the ellipsis . . . stands for operators that are in- 
serted at the beginning and the ends of paths and cor- 
respond physically to density, current, etc. Lt is the 
number of A nodes (or B nodes) along the vertical di- 
rection, interpreted as discrete (imaginary) time. (The 
operator (UbUa) t was denoted by U in Ref. 0). The 
supertrace STr implements periodic boundary condition 
along the vertical direction. Let there be Lx nodes in 
each row (the 'horizontal', or X direction). The opera- 
tor Ua (rcsp. Ub) is formed by multiplying all the Lx 
transfer matrices at A (resp. B) nodes in a given row. 
These matrices act in a tensor product of bosonic and 
fermionic Fock spaces defined for each vertical column 
of links. The quantum state in the Fock spaces in each 
horizontal row of links should be thought of as resulting 
from the quantum state of the previous (in discrete imag- 
inary time) row upon action of either Ua or Ub, within 
a single time step. The presence of a fermion or boson 
on 9-link represents an element of a path traversing that 



Both bosons and fermions arc needed to ensure 



the cancelation of contributions from closed loops (this 
ensures that STr (U b U a ) Lt = 1). 

Usually one needs two types of bosons and fermions, 
'retarded' and 'advanced', to be able to obtain two- 
particle properties (that is averages of products of re- 
tarded and advanced Green's functions) relevant for the 
calculation of transport properties. However, the sym- 
metry relation in Eq. ( |l6[) relates retarded and advanced 
Green's functions. Hence, as it turns out, for the compu- 
tation of low enough moments of, say, retarded Green's 
functions, we need only one fermion and one boson per 
spin direction per link of a row. (This will be discussed 
in detail below.) Wc denote them by fa, b a for the 
links going up (up links), and f a , b a for the down links. 
On the up links, f a , b a satisfy canonical fermion and 
boson commutation relations, respectively. But to en- 
sure the cancellation of closed loops, we must either take 
the fermions on the down links to satisfy the modified 
anti-commutation relations: {fa, fa} = ~3a,a' (while the 
bosons on the down links, b a and w a satisfy the canoni- 
cal commutation relations), or, alternatively, we take the 
bosons on the downlinks to satisfy the modified commu- 
tation relations: \ba,b\,\ = —5 a ,a' (while the fermions 
on the down links, f a and /J satisfy the canonical anti- 
commutation relations). Each node transfer matrix can 
be written in terms of these bosons and fermions. For any 
given disorder realization, this node transfer matrix turns 
out to commute with all the generators of the sl(2|l) 
superalgcbra, as shown in Rcf. [l7] (see Appendix A for 
a summary of the relevant superalgebra representation 
theory). Performing the average over disorder (indepen- 
dently on each link) projectsEH the Fock space of bosons 
and fermions into the fundamental (dual-fundamental) 3 
dimensional representation of sl(2|l) on up links (down 
links). These are precisely the states which are singlets 
under the random ('gauge') SU(2) on the links. 



D. Network boundary 

The extension of the above formalism to the case of 
a network with a reflecting boundary in the vertical di- 
rection is straightforward. Relegating some technical de- 
tails to Appendix [B, the reflecting boundary is seen to 
preserve the full sl(2 1) SUSY present in the bulk. Using 
the fully intact supersymmetry one can then retain the 
mapping to the perimeters of percolation clusters (we re- 
fer to them as hulls henceforth) obtained in Ref. D], but 
now the percolation hulls are confined to the half plane 
in 2D with a reflecting boundary (sec Fig. [Q). 

We will consider below correlation functions of oper- 
ators in the network model in the vicinity of the SQH 
transition.. When using continuum notation one needs 
to recallc3Sa that, in general, operators acquire addi- 
tional singularities when approaching a boundary; they 
may vanish or diverge. By saying that points lie on the 
boundary, we imply that these points are not literally lo- 
cated at the boundary, but rather that the distance of the 
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points from the boundary is much smaller than the dis- 
tance of these points from each other, and much smaller 
than the system size. The same general reasoning will 
apply when we consider multifractality of wave functions 
at the boundary. Let us recall that in the literature of 
surface transitions in ordinary magnets, various kinds of 
cases which are often referred to as 'ordinary, 'extraor- 
dinary' and 'special' boundary transitons are discussed, 
referring to the respective ordering of bulk and bound- 
ary. The case of the reflecting boundary condition that 
we consider here is analogous to the so-called 'ordinary' 
surface transition where the boundary and the bulk un- 
dergo the LD transition simultaneously. 

We will need the bulk and boundary scaling dimen- 
sions of the 1-hull operator of critical percolation in var- 
ious calculations below. This has ber.p derived using a 
variety of techniques in the litcratureHO'EHl and is found 
to be x s = 1/3 for the boundary 1-hull operator. The 
corresponding bulk field has scaling dimension xi, = 1/4. 



IV. MULTIFRACTALITY USING 
SUPERSYMMETRY 

A. Wave function correlators 

In the context of LD transitions, multifractality man- 
ifests itself through the anomalous scaling with system 
size of the moments of the square of the wave function 
amplitude, |^>(r)| (wave function ' intensity') at critical- 
ity (see, for example, Refs. ||and 20). The wave function 
moments conventionally appear in the form of the so- 
called averaged generalized inverse participation ratios 
(IPR) P q (the overbar denotes disorder average) 



P„ 



Ma 



d u *r x |^(r x )| 2 ?ocL^|V(r K )| 2 ?, (17) 



whose scaling with linear system size L at a critical point 
defines the critical exponents r q and (often referred 
to as 'anomalous dimension'), 



P q oc L- T - 



= dq 



A* 



qHx ~ D x , 



(18) 



Here, D x is the dimension of the region M x over which 
the integration is performed in (jl7j). In particular, we 
will be interested in the cases of bulk (x = b), boundary 
(x = s) and corner (x = c) where D x takes on, respec- 
tively, the values d, d — 1 and d — 2. The corresponding 
anomalous dimensions will be denoted by A* A^ and A£. 
One notes that A\ = due to the fact that wave func- 
tions are normalized to unity with respect to the whole 



system (implying P q= \ — 1 and (if, = 0). There is no 
such general constraint for A\ and AJ. However, we 
adopt the convention of setting them equal to zero and 
the corresponding contribution to rf is reflected in fi s 
and /i c respectively. With this choice, a nonvanishing ex- 
ponent fi x characterizes a non-trivial scaling behavior of 
the local wave function intensity at the sample location 



r x , 



L d |^(r x )| 2 oc 1/ZA 



(19) 



which is known to occur in certain iio&i-convcntional 
Anderson localization symmetry classesBH3, as was dis- 
cussed previously in Ref. ||(| 

The wave function moments have a natural interpre- 
tation as operators in a field theory describing the LD 
transition. As an aside, it is worthwhile emphasizing 
that the interpretation of multifractal moments in terms 
of field theory operators is a very general one and not just 
restricted to LD transitions (see Ref. |22| for a discussion 
of this point). We have defined multifractal exponents 
in terms of moments of wave function at a single point 
in the sample in Eq. (17). In field theory, for calculat- 
ing the multifractal exponent A^, (restricting q to be an 
integer,) it is easier to find the scaling of correlation func- 
tion of q operators, each corresponding to the wave func- 
tion intensity L 2 \i/j(r)\ 2 . The g-point correlations of crit- 
ical eigenfunction intensities take on a convenient scaling 
form when the distances between all points are equal to 
the same scale r, that is |r,- — r,|~ r and it follows that 



L^+^)|V>(ri)V(r 2 )...VM 5 



(20) 



where L is the linear system size. 

Off criticality (that is, for e / or 7 > in the 
SQH system), multifractal scaling behavior holds oh. 
length scales much shorter than the localization lengthEZl 
£ 7 ~ 7~ 4 / 7 , beyond which the wave function amplitudes 
are exponentially smal£jc3. This implies that off criti- 
cality, the above q-point correlator should be written as 



(21) 



Assuming that all wave functions at a given energy show 
statistically identical multifractal scaling behavior, one 
can write the L.H.S of the above correlator as 



■ ip(r q )\ 2 



En.^.i, IV'n^i)^^) ■ ..ipr q {r q )\ 2 8{<n - e ll )S(e 2 -e t2 ).--S(e q - a J 



p(ei) p(e 2 ) • • • p(e q ) 



(22) 
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where i\, i 2 . . . i q are additional quantum numbers required, besides position r, to label the wave function (energy and 
spin in the SQH case), p represents the global density of states. Taking ei, e 2 . . .e q ~ e and jr, — rj\ ~ r , we see that 
the exponent is given by the scaling behavior of V q (ri,r 2 . . .r g ;e) with respect to r, where the function V* is 
defined by 



T> x q ( ri> r 2 ...r g ;e)= ]T |^ (n)^ (r 2 ) . . . V% (^)| 2 5(e - e^Sfe - e la ) . . . 6(e - e iq ). (23) 

21 ,12- ..iq 



We understand that the points ri . . . r q are all chosen to 
lie in the region M. x considered in Eq. (p"T 

we 



Using the Green's functions defined in Sec. [II B 



can write the above wave function correlator (ignoring 
overall factors of (27r)~ 9 ) as 



i ) q{r 1 ,r 2 . . -r q ; z) oc tr [G R (r k ,r k ) - G A (r k ,r k ) 



k=l 



(24) 



Here 'tr' denotes the trace over the spin indices. As was 
discussed in Rcf. the multifractal exponents can also 
be obtained from a different Green's function product. 



V q (n,r 2 ...r q ;z)(x tr JJ[Gij(r fe , r k +i) - G A (r kl r k+1 )}, 
fc=i 

(25) 

where r q+ \ = r%. When |r,; — Tj\ ~ r , this product shows 



the same scaling behavior, r « as the previous one. Since 
the calculations of both functions V and V are very sim- 
ilar, we have explained in detail only the calculation of 
V. By suitably choosing the correct SU(2) invariant (see 
next section), we will extract below the scaling of V. Our 
results agree with those previously obtained in Ref. [50| for 
the bulk case, x = b, but we obtain these results using a 
different method (namely that developed in Ref. [ll], using 
supersymmctry). For the other s (boundary) 

and x = c (corner), our results are entirely new. 



B. Calculation of A? 



We start with the calculation of the exponent for 
the case q — 2. The corresponding Green's function 
product that we have to evaluate is (dependence on z = 
= exp(— 7), defined in Eqs (|l~4| , |l~5| ), is understood) 



VI (n , r 2 ) cx tr [G R (n , n) - G A (n , r x )\ tr [G R (r 2 ,r 2 )-G A (r 2 , r 2 )] . (26) 
Using Eq. @, we can write the above relation in terms of retarded functions alone as 

VI (n, r 2 ) oc tr [G R (n, r x ) - 1] tr [G R (r 2 , r 2 ) - 1] 

= trGfl(ri,n)trG fl (r 2 ,r 2 ) - tr G R {r l7 n) - tr G fi (r 2 , r 2 ) + 1. (27) 



We will use the supersymmetry teehniqueu to implement 
the disorder average and expressEZI all averaged Green's 
functions using the second quantized formalism in terms 
of generators of the sl(2|l) Lie superalgebra. (A sum- 
mary of certain basic elements of the relevant represen- 
tation theory of the sl(2|l) superalgebra is presented in 
Appendix [A|.) 

The retarded Green's functions can be expressed as 
expectation values of canonical boson and fermion oper- 
ators. (From here on, all Green's functions are under- 
stood as retarded unless they have an explicit subscript 
^4). Specifically, we can write the following Green's func- 
tions, in any fixed realization of disorder, as (the sub- 
scripts below refer to both, the spatial coordinate and 



the spin index) 

G tj = STr[b i b J ((U B UA) LT ]=(b i b]), 
G kn = STr[f k P n (U B U A ) LT ] = (fjl), (28) 
G l3 G kn = STr[b t b}f k fl(U B U A ) L -] = (b^fjl). (29) 

(Here, the fact that STr (UbU a ) Lt = 1 was used.) Let us 
take i -> (ri,ai), j -> (n,0i),k -> (r 2 ,a 2 ),n -> (r 2 ,/3 2 ) 
in the above equations. The LHS of the last equation thus 
reads G^ 1 (ri,ri)G^ (r 2 , r 2 ). After performing the disor- 
der average, only SU(2) singlet combinations are non- 
vanishing. Hence we contract both sides with the SU(2) 
invariant S^S^ giving LHS = tr G(n, ri) tr G(r 2 , r 2 ). 
The RHS of Eq. (^9|) is (we can raise and lower indices 



8 



trivially since the metric is the unit matrix) 

RHS = S%6% (Mn)^>i)/a 2 (r 2 )4 2 (r 2 )) 

= ((2B(ri) + l)(l-2Q 3 (r 3 ))>. (30) 

(The expressions for the generators B and Q 3 of the Lie 
superalgebra, as reviewed in Appendix A, were used.) 
Thus the LHS and RHS together imply: 

trG(n,ri)trG(r 2 ,r 2 ) 

= ((2B(ri) + l)(l-2Q 3 (ra))). (31) 

Similarly, considering Eq. (p8|), let k — > (r 2 ,a 2 ),n — > 
( r 2,/3 2 ) giving LHS = (r 2 ,r 2 ). To do disorder av- 
erage, contract with the SU(2) invariant 8%l, giving 
LHS = tr G(r 2 ,r 2 ). The RHS of Eq. © gives 

RHS = 5% (/ Q2 (r 2 )4 2 (r 2 )) = <1 - 2Q 3 (r 2 )}. (32) 
Hence from LHS and RHS, we obtain 

trG(r 2 ,r 2 ) = (l-2Q 3 (r 2 )). (33) 

Now using Eq. @ and Eq. @, we can write Eq. @ 

as 

P 2 K (r 1 ,r 2 )«-(B(r 1 )Q 3 (r 2 )). (34) 

We will now compute this correlator off criticality, at a 
finite correlation length £ 7 , arising from a non- vanishing 
'broadening' 7 > 0, that is z = cxp(— 7) < 1. 

The angular brackets in Eq. (|34|) denote the supertrace 
STr taken over the full Fock space of canonical bosons 
and fermions on each link of the network before disorder 
averaging. As explained in Ref. [l?], the disorder aver- 
age projects this Fock space into the fundamental three- 
dimensional representation of sl(2|l) on the up links and 
the corresponding dual representation on the down links. 
To sum only over the three states in this representation, 
we use the notation 'str' (see Eq. (|35|)). We note that 
the state with odd fermion number iu both of these rep- 
resentations has has negative normtll. The prefix 's' in 
str denotes the fact that these negative norm states con- 
tribute to the trace with a negative sign. This is es- 
sential to get the correct overall factor of unity for each 
loop traverscdO. The action of the node transfer ma- 
trix, Tsu, on the tensor product of the disorder-averaged 
states on an up link and a neighboring down link can 
be represented by a linear combination of the only two 



sl(2|l) invariant operators in this 3x3 dimensional vector 
space. These are the identity operator and the projec- 
tion operator onto the singlet state, with weights 1 — t 2 Sa 
and t 2 Sa respectively. These weights can be considered as 
the classical probability of a state turning left or right 
at a given node under the action of the node transfer 
matrix. When we multiply all such node transfer matri- 
ces together, the partition function can be represented as 
a sum over all classical configurations of densely packed 
loops (the closed classical paths along which the disor- 
der averaged states propagate) on a square lattice. Each 
loop gets an overall weight which is the product of the 
probabilities of turning in a particular direction at each 
node. The loops can be interpreted (sec Ref. [Tt] for a 
diagrammatic perspective) as the external perimeters of 
a cluster percolating along the bonds of a square lattice. 
This completes the mapping to percolation identified in 
Ref. 0. 

Now, in order to evaluate a disorder averaged correla- 
tor such as the one in Eq. ([34j), consider the 3 dimen- 
sional representation on the links (for example, the case 
where n, r 2 lie on up links) and all calculations are done 
in this representation. A loop passing through a link 
corresponds to all three states in the three dimensional 
sl(2|l) representation propagating on that link. Away 
from criticality, we assign a factor of z 2 ( B +Q 3 ) for each 
such link since the operator 2(B + Q 3 ) counts the number 
of states propagating on that link. The product of these 
factors along a path on the network gives the same power 
of z that occurs in the Taylor expansion of a matrix ele- 
ment of the Green's function in Eq. (jl4|). After multiply- 
ing such factors for all links through which a given loop 
passes, we also multiply it with the overall weight coming 
from the classical probability of the loop turning at each 
node as mentioned in the previous paragraph. Further, if 
there are operators inserted at specific points on the lat- 
tice, the only loops contributing to the partition sum are 
those which are constrained to pass through these points. 
Considering the correlator in Eq. ( |34"[ ) , there are two dif- 
ferent kinds of loop configurations which contribute to 
the partition sum: (1) a single loop passes through both 
points ri and r 2 , (2) two different loops pass through 
each of these two points, T\ and r 2 . These two terms are 
the probabilistic versions of the usual connected and dis- 
connected parts of any correlation function. Writing the 
contribution for each of these types separately and sum- 
ming over all possible loop configurations together with 
their respective weights, we write Eq. ( p4| ) as 



4(B(n)Q 3 (r 2 )) 

(1 o\ A \ /-i o\ A \ 

= V str 2 z 2N ^ z 2N ^ \ P(n, r 2 ; N 12 , N 21 ) 

n 12 ,n 21 \0 1/ \0 z 2N ^ J \ I \Q z 2N ^ ] 
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(i o o\ (i o o \ /— i o o\ A o o \ 

+ V str 2 z 2N str z 2N ' P-(n, r 2 ; AT, N') 

n,n> \0 1/ \0 z 2Ar / \ I J \o z 2JV 7 

= " E C 1 - P(n,rr,N) £ [l - z 2W ] [l - z 2W '] P_ (n ,r 2 ;N, N'), (35) 

JV 2V,iV' 



where P(ri,r 2 ; N) is the probability of a single loop 
of length ./V passing through both r\ and r 2 , and 
P_(n, r 2 ; AT, AT') is the probability that in a given perco- 
lation configuration the points r\ and r 2 belong to two 
different, non-overlapping loops of lengths A" and N', re- 
spectively. In the above equation, the probability fac- 
tor P(ri,r2\ N12, N21) is for a path of length N 12 going 
from r± to r 2 and then a path of length N 2 i returning 
from r 2 to n. In the next line, we simplified this using 
the fact that this is the same as a loop of overall length, 
A" = N12+N21 passing through both the points, and set- 
ting P{r 1 ,r 2 \ N) := P(ri,r 2 ; N 12 , N 2 i)- The appearance 
of these classical percolation probabilities in the above 
formula arises from the product of the individual fac- 
tors for turning left or right at each node of the network 
model. Thus, using Eq. ([35]), we can write Eq. (|34|) as: 



N s/7 N i-8/7 



(38) 



V*(n,r 2 ) <x^[l- z 2N ]P(r u r 2 ;N) 

N 

+ E [l-z 2N }[l-z 2N ']P-(ri,r 2 ;N,N'). 

(36) 



N.N' 



We see that the leading order terms turn out to van- 
ish in the critical limit z — > 1, as observed in Rcf. |3^. 
So we need to consider the sub-leading behavior. The 
percolation probabilities discussed above are the Laplace 
transformed versions of the correlation functions of 1-hull 
operators in percolation. The 1-hull operator represents 
the density-of-states operator in the SQH problemtH. Wc 
know the scaling dimension of these o perat ors as men- 
tioned in the last paragraph in Sec. [IIP . Under the 
Laplace transform, the variable N is the conjugate of the 
energy e + ij = — i In z which is represented by the 1- 
hull operator. Now, using the scaling dimensions of the 
1-hull operator to deduce the scaling of the Laplace trans- 
formed correlation functions (see, for example, Ref. 51), 
one finds the following leading scaling behavior of the 
above-mentioned probabilities 



P(n,r 2 ;N) ~ N- 2/(2 - Xb) r- Xb 
~ TV-s/V- 1 / 4 , 
P_(r 1; r 2 ; N, N') ~ P(n; N)P(r 2 ; N') 



(37) 



where P(n;N) is the probability that the point n be- 
longs to a loop of length N. The points r\ and r 2 lie 
in the bulk and r = \n - r 2 \ < £ 7 ~ 7~ 4/7 ~ N 4 / 7 . 
Similarly, the scaling expressions in the case where both 
points lie on the boundary are 

P(ri,r 2 ;AT) ~ N -i-x s /(2-x b ) r -x a 

„ JV -25/21 r -l/3 j (39) 

P-Cri.ra; N, N') ~ P(n; N)P{r 2 ; N') 

^ JV -2S/21j V -/-25/2X _ ( 40 ) 



We see that only the first term in Eq. (|36j) (which is 
the connected part) gives rise to the non-analytic r de- 
pendence in the limit r <C £ 7 , which we are interested in 
computing. Hence the scaling of T> 2 {r\,r 2 ) is given by: 



N 

fr" 1 / 4 (bulk) , 

I r — 1/3 (boundary). 



(41) 



Thus, upon comparison with Eq. (|2lD, we finally read off 
the multifractal exponent for q = 2 as being = —1/4 
for the bulk, and A 2 = —1/3 for the boundary. We will 
discuss the case of the wave function scaling behavior at 
corners at the end (see Sec. VI) as it is a straightfor- 
ward extension of the boundary case upon making use of 
conformal invariancc. In the following paragraph we con- 
sider the more interesting exponent describing the staling 
of the third moment of the square of the wave functions 
amplitude. 



C. Calculation of A3 

The algebraic procedure for calculating the multifrac- 
tal exponent A 9 for q = 3 is almost same as that we used 
above for q = 2. Hence we present only the important 
steps and focus on the main results and the interesting 
differences. We start with the expression: 



f> x (ri,r 2 ,r 3 ) cx tr [G R (ri,n) - GU(rL,n)]tr [G R (r2,r 2 ) - Ga^, r 2 )] tr [G R (r 3 ,r 3 ) - Ga(t 3 , r 3 )]. 



(42) 
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Converting all advanced Green's functions to retarded ones using Eq. (|16|), we can write the above equation as: 

V%{n,r 2 ,r 3 ) cx tr[G(n,ri) - 1] tr [G(r 2l r 2 ) - 1] tr [G(r 3 ,r 3 ) - 1] 

= trG(r uri )trG(r 2l r 2 ) trG(r 3 ,r 3 ) - 1 - ]T trG^, n) trG(r 2 , r 2 ) + ^ trG^n), (43) 

[ri,r 2 ,r 3 ] [ru^rs] 
I 



where Y~\ i denotes sum over terms with all cyclic 
permutations of (ri, r 2 , r 3 ). 

The new piece we have to evaluate here is the product 
of three Green's functions. For this consider the following 
two identities, obtained by applying Wick's theorem, and 
valid in any fixed realization of disorder: 

[G^Gim + G im Gij] G kn = (hb]b t blf k fl) . (44) 



Exchanging bosons and fermions in the above equation, 
we find 

[GijG lm - Gi m Gij] G kn = (fif]fifib k bl). (45) 

In Eq. @), let i -> (ri,ai), j -> (ri,/3i), I -> (r 2 ,a 2 ), 
m -> (r 2 ,f3 2 ), k -> (r 3 ,a 3 ), n -> (r 3 ,/? 3 ), yielding 



LHS = [G;;(r 1 ,r 1 )G; 2 2 (r 2 ,r 2 ) + G; 2 1 (ri,r2)G;;(r 2 ,r 1 )]G; 3 3 (r 3 ,r 3 ). (46) 
As above, it is convenient to contract with the SU(2) invariant , after performing the disorder average 



tr G(n, ri) tr G(r 2 , r 2 ) tr G(r 3 , r 3 ) + tr G(r 3 , r 3 ) tr [G(n, ra)G(ra,ri)] 



(6a(ri)6i(ri)6^(r a )6j(r a )/ 7 (r 3 )/t(r 3 )) = ( [1 + 2B( n )} [1 + 2B(ra)] [1 - 2Q 3 (r 3 )] ). 



(47) 



We see that the first term on the LHS of Eq. (|47|) is the Green's function product we want to evaluate in Eq. (|43]), 
while the second term has to be eliminated. This can be achieved by performing steps, analogous to those above, in 
Eq. @, yielding 



tr G( ri , ri ) tr G(r 2 , r 2 ) tr G(r 3 , r 3 ) - tr G(r 3 , r 3 ) tr [G(r u r 2 )G(r 2 , ri )} 
= </«(ri)/l(r 1 )/ /3 (r 2 )4(r 2 )6 7 (r 3 )6t ( r3) ) = ( [i _ 2 Q 3 (r 1 )] [1 - 2Q 3 (r 2 )] [1 + 2B(r 3 )] >. 



Adding Eq. ([47]) and Eq. we obtain 



2trG(r 1 ,r 1 )trG(r 2 ,r 2 )trG(r 3 ,r 3 ) = ( [1 + 2B(n)] [1 + 2fl(r 3 )] [1 - 2Q 3 (r 3 )] 

+ [1 - 2Q 3 (n)] [1 - 2Q 3 (r 2 )] [1 + 2B(r 3 )] ). 

We can now use Eqs. ||[ |49|) to write Eq. (|43j) as 

^(n,r 2 ,r 3 ) cx ( [1 + 2S(n)] [1 + 2B(r 2 )} [1 - 2Q 3 (r 3 )] + [1 - 2Q 3 (n)] [1 - 2Q 3 (r 2 )] [1 + 2S(r 3 )] ) 
- <[l-2Qs(ri)][l + 2B(r 3 )]>+ £ <1 + 2B(n)> - 1. 



[ri,r2,r 3 ] 



[fi,r2,r- 3 ] 



(48) 



(49) 



(50) 



Converting these expressions to percolation probabilities 
is exactly analogous to the two point case. Taking the 
critical limit z — ► 1, we obtain 



given by that of the usual 3-point correlation function of 
percolation 1-hull operators at criticality, 



i>U r i> r 2> r 3) P(ri,r 2 ,r 3 ) + P(r 1 ,r 3 ,r 2 ). (51) P(n, r 2 , r 3 ), P(r x , r 3 , r 2 ) 



-3a:i, — 3/4 



3 x s _ , 1 



where P{r\, r 2 , r 3 ) is the probability of a loop of any size 
traversing ri,r 2 ,'r 3 in that order. There is no cancela- 
tion to leading order at criticality here, unlike the two 
point case. The final scaling of this correlator is simply 



(bulk), 
(surface). 



(52) 



Hence the value of multifractal exponent for q = 3 is 
A 3 = —3/4 and A| = —1, for bulk and surface, respec- 
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lively. 



D. Higher multifractal exponents 



The procedure for calculating A| was very similar to 
that of Af (albeit more tedious) and one might ask if 
it can be extended to higher multifractal exponents, 
with q > 3. On the other hand however, the fact that 
we are able to calculate A 3 at all using our formalism is 
conceptually surprising because usually in supersymmct- 
ric problems, an additional set (='replica') of boson and 
fermion operators is required for each additional Green's 
function factor entering the product in Eq. (p4|). Here we 
are able to calculate both two and three-point functions 
with the same number of boson and fermion operators 
(replicas) . 

To understand this, we look at the calculation of A 3 
carefully. We had an extra 'unwanted' Green's function 
product in Eq. (f47|). But we were able to eliminate it 
by using an equation similar to Eq. (|48|), with bosons 
and fermions exchanged. The unwanted Green's func- 



tion product canceled between Eq. (47) and Eq. ( |4g| ) 
when added together, giving us the exact product that 
we wanted. 

Now does such a cancelation go through for higher 
point functions? To answer this, we look at the next 
higher exponent, Af . Here we will have to evaluate a 
product of four Green's functions, 

trG(ri,ri)trG(r 2 ,r 2 )trG(r 3 ,r 3 )trG(r 4 ,r 4 ). 

To evaluate this, we will have to use the following identity 
obtained from Wick's theorem (in a fixed realization of 
disorder): 



[GijGlm GimGlj] \GpqG rs ~t~ Gp S G r q\ 

= {fif]fiflh P b\b r b\). 



(53) 



In this equation, let again i — > (r l7 ai), j — > (ri , /3i), 
/ -> (r 2 ,a 2 ), m -> (r 2 ,/? 2 ), P -> (r3,a 3 ), q -> {r3,0a), 
r — ► {ri, a.4), s — > (r 4 ,/? 4 ). Then contracting with the 
SU(2) invariant 5&<5****& 

gives us upon averaging the 

relation: 



tr G(n, n) tr G(r 2 , r 2 ) tr G(r 3 , r 3 ) tr G(r 4 , r 4 ) + tr G(n, n) tr G(r 2 , r 2 ) tr [G(r 3 , r 4 )G(r 4 , r 3 )] 
- tr [G{r u r 2 )G{r 2 ,r 1 )} tr G(r 3 , r 3 ) tr G(r 4 , r 4 ) - tr [G{r u r 2 )G{r 2 ,n)\ tr [G{r 3 , r4: )G(r 4 ,r 3 )} 

= ( [1 - 2Q 3 (n)] [1 - 2Q 3 (r 2 )] [1 + 2B(r 3 )} [1 + 2B(r 4 )] ). (54) 

There are three unwanted terms on the LHS. Exchanging bosons and fermions gives on the other hand: 

tr G(n , n) tr G(r 2 , r 2 ) tr G(r 3 , r 3 ) tr G(r 4 , r 4 ) - tr G(n , n) tr G(r 2 , r 2 ) tr [G(r 3 , r 4 )G(r 4 , r 3 )] 
+ tr r 2 )G(r 2 , n)] tr G(r 3 , r 3 ) tr G(r 4 , r 4 ) - tr [G(ri, r 2 )G(r 2 , n)] tr [G(r 3 , r 4 )G(r 4 , r 3 )] 

= ( [1 - 2Q 3 (r 3 )] [1 - 2Q 3 (r 4 )] [1 + 2B(n)] [1 + 2B(r 2 )] ). (55) 

Adding Eq. (|5^) and Eq. (^) eliminates two unwanted pieces: 

2 tr G(r 1? n) tr G(r 2 , r 2 ) tr G(r 3 , r 3 ) tr G(r 4 , r 4 ) - 2 tr [G(ri, r 2 )G(r 2 , n)] tr [G{r 3l r 4 )G( ri ,r 3 )] 
= ( [1 + 2fi(n)] [1 + 2B(r 2 )] [1 - 2Q 3 (r 3 )] [1 - 2Q 3 (r 4 )] ) + 
( [1 - 2Q 3 (n)] [1 - 2Q 3 (r 2 )] [1 + 2B(r 3 )\ [1 + 2B(r 4 )] >. (56) 

I 



We are still left with one unwanted piece which cannot 
be evaluated or canceled with something else. Any other 
combination of the generators of the Lie superalgebra at 
four points will generate more terms on applying Wick's 
theorem and cannot be canceled out. 

Now we see the special feature of the three point cal- 
culation. We had two Green's function products (by ap- 
plying Wick's theorem), of which only one was necessary. 
The supersymmetric formulation gave us one more equa- 
tion due to boson-fcrmion intcrchangcability. The un- 
wanted piece canceled between the fcrmionic and bosonic 
equations. This fact does not help us in higher correla- 



tion functions as they have more unwanted pieces. It is 
also clear that situation becomes worse for higher n-point 
functions. Interestingly, the same conclusion was reached 
in a very different way in Ref. |30|. 



V. LOCAL DENSITY OF STATES AND POINT 
CONTACT CONDUCTANCE 

Having considered the multifractal calculation in de- 
tail, the calculation of other boundary critical exponents 
is completely analogous and we list only the important 
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results. Some of the bulk exponents have been found in 
Refs. HI and H 

using a very different technique. 
The averaged local density of states (LDOS), summed 
over the spin indices, can be written in terms of Green's 
functions as 



1 



(p*(r,e)) = — tr [G R (r,r) - G A (r,r)} 



— trG fl (r,r)-l. 



(57) 



This can again be expressed in terms of the sl(2|l) super- 
symmetry generators as (l/27r)^2B(r)), and thi s ave rage, 



following the same steps presented in Section IV B 



be written in terms of percolation probabilities as 
(2B(r)) =1~Y1 P ( r > N ) cos 2Ne - 

N 



(58) 



(The same result was also obtained in Rcf. 46 using, as 
mentioned, different techniques.) As we have mentioned 
in the end of Sec. Ill, the boundary and the bulk scaling 
dimensions of the one hull operator are x s = 1/3 and 
xt, = 1/4. The latter value implies that the percolation 
hull has fractal dimension 2— Xf, = 7/4, so that P(r, N) ~ 
f or r m bulk. This yields, according— to Eq. (57), 
the following scaling behavior of the LDOSLJ 



p\e) cx e ^/(2-^) = e 1 ' 7 . 



(59) 



Note that 2 — x^ = 7/4 is the dynamic critical exponent 
governing the scaling of energy with the system size L 
at SQH criticality, so that the level spacing at e = 
(and thus the characteristic energy of critical states) is 
6 ~ L- 7 / 4 . In our case, when the point r is located at 
the boundary we find 



P(r,N) ~ A" 



-l-x a /(2-x h ) 



N' 



-25/21 



and the LDOS scaling 

p s (r, e )(X£ I '/ (2 -^= e 4 / 21 . 



(60) 



(61) 



Here we have used the bulk dynamic critical exponent in 
determining the energy scaling of the boundary LDOS. 
This is because we are dealing with an 'ordinary sur- 
face transition' and here surface (boundary) criticalitjcJs 
driven by the divergence of the bulk correlation lengthEil. 
Thus the scaling of LDOS changes between bulk and 
boun dary. This means, as was mentioned in Section 
EVA, that the average of the (square of the) wave func- 
tion amplitude is suppressed at the boundary, giving rise 
to a non- vanishing value of fi x=s = 1/3 — 1/4 = 1/12 (see 
Eq. ©). 

A similar procedure is adopted for the calculation of 
other boundary exponents (we give a table of all expo- 
nents later). The boundary diffusion propagator can be 
written as 



(n M (n,r2)) = 2(V_(ri)W + (r2)>, 



where r\ and r 2 lie at the boundary. In terms of perco- 
lation probabilities this reads 



{n ss (r 1 ,r 2 ))=2^2P(r 1 ,r 2 ;N)z 



2.V 



N 



Taking the limit z — > 1 (critical point), gives 
<n ss (n,r 2 )) =2j2P(n,r 2 ;N) = 2P(n,r 2 ] 

7V7 



(63) 



(64) 



Here the probability P(ri,r 2 ) for the points riand r 2 
to be connected by a hull of any length scales at the 
boundary as r~ 2Xs 



giving 



(ILj a (ri,r 2 )) ~ |n - r 2 \ 



-2/3 



(65) 



(62) 



Another physical quantity of interest is the boundary 
point-contact conductance. In a network model setting, 
this is defined as the conductance between two boundary 
links n and r 2 which are cut to make it possible to insert 
and extract currents from thcrnEE In the second quan- 
tized formalism, this is equivalent to creating a boson (or 
alternately a fermion because of the supersymmetry) for 
each spin direction at one link and destroying the same 
particle at another link. This translates to the expression 

(9 PO int(ri,r 2 )) ee (fiinWiin)^)^)}. (66) 

Iii terms of sl(2 1 1 ) generators, this equals 
(Q+(ri)Q-(r 2 )). Skipping here the mapping to 
percolation probabilities, we find that the point-contact 
conductance decays exactly in the same way as the diffu- 
sion propagator, and, in particular, scales as \r\ — r 2 \~ 2 ^ 3 
at the boundary, at criticality. This result is expected 
since the correlators for the diffusion propagator and 
the point-contact conductance represent different su- 
pcrspin components in the same representation of the 
superalgcbra sl(2|l). 



VI. EXPONENTS IN OTHER GEOMETRIES 

There are two distinct ways of extending our discussion 
of boundary behavior. The first one is when some of the 
points lie on the boundary while the others lie in the 
bulk. The other case is to consider boundaries with more 
complicated geometries, the simplest example of which 
would be a wedge with opening angle 9. (The boundary 
case, considered in the previous sections corresponds to 
9 = ?r). 

The two-point quantities are easily computed when one 
point T\ lies in the bulk and the other point, r 2 is at 
the boundary. These scale as r~( ?,b+ ' , ii)/ 2 where r\b and 
r?|| are the usual exponents giving the decay of the two 
point function in the bulk and along the boundary, re- 
spectively. Hence the diffusion propagator and the point- 
contact conductance between a point in the bulk and 
another at the boundary, both scale with distance r as 

r -7/12_ 
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In the case of multifractal exponents, we can calculate 
the scaling behavior of correlation functions similar to 
those in Eq. (|24|). But we should not interpret these as 
representing properties of a single multifractal since mul- 
tifractality is essentially a single point property. With 
this caveat, the value of the quantity analogous to A2 
when one point is in the bulk and another is at the bound- 
ary is —1/4. 



Geometry 


O 


• 




LDOS 




£ 4/21 


e (4/21)(^/ 9) 


A 2 


-1/4 


-1/3 


-7r/3(9 


A 3 


-3/4 


-1 


-7r/<9 



TABLE I: One-point exponents in various geometries: bulk, 
near straight boundaries and near corners. The first line 
indicates the scaling behavior of the local density of states 
(LDOS) at the SQH transition. The second and third lines 
represent the multifractal scaling exponents for the second 
and third moments of the critical wave function intensity. 



Geometry 










f 


II(ri, r 2 ) 


r -V2 


r -7/12 


r -2/3 


r -l/4-(l/3)(vr/e) 



TABLE II: Critical scaling behavior of the diffusion propaga- 
tor, n(n,r2) in various geometries. The point-contact con- 
ductance has an identical scaling behavior as indicated in the 
text. 

The next generalization is to analyze the SQH transi- 
tion in a wedge geometry. This can be .readily carried out 
using conformal invariance arguments^. The conformal 
transformation w = z 6 ^ maps the boundary geometry 
in the z-plane to a wedge of opening angle 9 in the w- 
plane. From this, it can be deduced that if we consider 
a 2-point function with one point lying near the wedge 
tip and another deep in the bulk at distance r, the two 
point function decays as r~ ne where 



V6 



7T 



(67) 



This enables us to calculate the relevant exponents. In 
the density of states calculation, we will have to replace 
x s by (ir/0)x s . 



VII. CONCLUSIONS AND OUTLOOK 

The central result of this paper is the calculation of 
boundary critical and multifractal exponents for the spin 
quantum Hall transition in various geometries (see Tables 
p| and p| for a summary of results). In this paper, we have 



confined ourselves to the exact calculation of anomalous 
multifractal dimensions for low moment order q, using 
the percolation mapping. In general, in the theory of 
multifractality, the information contained in the set of 
all the critical exponents defined in Eq. ( |l8| ) can also 
be expressed in terms of the so-called singularity spec- 
trum f x {oi), its Legendre transform. We remark that 
in the presence of a boundary, the singularity spectrum 
concept has to be understood more broadly. This is be- 
cause unlike in ordinary critical phenomena, the presence 
of a boundary can affect the singularity spectrum (and 
t£) of the entire system, including bulk and boundaries, 
in a significant way, even in the thermodynamic limit. 
These issues were first pointed out in Ref. P6| and were 
analyzed therein, as well as in Refs. ^ and |53|. 

The analysis of multifractality in various other geome- 
tries leads to interesting new concepts. Some of these 
ideas haveJpeen explored for the Anderson transition in 
d = 2 + (£3, power-law random baaded matriceso and 
the 2D symplectic class transitionE§E3 all of which dis- 
play the characteristic multifractal property of LD criti- 
cal points. A, similar study of Dirac fcrmions in random 
gauge fieldsEa and the integer quantum Hall transitions 
in two dimensions is expected to further our understand- 
ing of boundary multifractality and provide important 
clues regarding the structure of some of the unknown 
bulk theories. 
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ary MF at the SQH transition and A. D. Mirlin for a 
previous collaboration on boundary MF. This work was 
supported in part by the Chandrasekhar and Sachs fel- 
lowships (ARS), NSF Career award DMR-0448820, NSF 
MRSEC DMR-0213745, the Alfred P. Sloan Foundation 
and the Research Corporation (IAG), and NSF grant 
DMR-0706140 (AWWL). 



APPENDIX A: REPRESENTATIONS OF sl(2 1) 
SUPERALGEBRA 



The sl(2|l) superalgebra has eight generators, of which 
four are bosonic (B, Q3, Q + , Q_) and four are fermionic 
(V+, V-, W+, W-). We use the convention of Ref. |6| for 
the generators. These satisfy the same commutation^ [, ]) 
and anticommutation ({, }) relations as the generators of 
osp(2|2) superalgebra: 



[B,Q 3 ] = [B,Q ± ]=0 
1 



[B,V ± ] = -V±, 
[Qa,Q±\ = ±Q±, 

4^ 



[Qh 



[Q3,V±] 

[Q+,V.] = V +[ 
[Q_,F + ] = F_ ; 
,V+] = [Q + ,W + 
{V+,V-} 



[B,W±] 
[Q+,Q-] 
[Qs,W±] 
[Q+,W-] 



1. 



[Q-. 
-- [Q-,V-] = 
{W+,W-} 



2 

2Q 3 , 

W+, 
W+] = W-, 
[Q-,W-] = 0, 

= 0, 
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{V + ,W+} = Q+, 
{V-,W-} = -Q-, 



{V+,W-} = B - Q 3 , 
{V-,W+} = -B-Q 3 . 



(Al) 



An important subalgebra is gl(l|l) formed by the gen- 
erators (B, Q 3 , V- , W + ). This is the SUSY present in the 
SQH network at finite broadening 7, when z = e -7 < 1 
(see Eqs. |l5|)), as well as at an absorbing boundary. 

The sl(2|l) algebra has an oscillator realization formed 
by bilinear combinations of the fermion and boson opera- 
tors on each link that are SU(2) singlets. For the up- links 
the oscillator representation is: 



Q 3 
Q+ 

B 
V+ 
V- 
W+ 
W- 



5(/ T V T + /!/;-!), 
l -(b\b ]+ b\b [ + l), 

^(Vi-W- 



(A2) 



These operators act irreducibly in the fundamental 
three-dimensional representation of sl(2|l) (denoted by 
3) with the space of states spanned by three SU(2) sin- 
glet states which we denote as \m), m = 0, 1, 2: 



vacuum) 



v + n = ±={b\f\- b \f\m, 



-|0> = /|/||0) 



(A3) 
(A4) 

(A5) 



We need the matrix elements of the sl(2|l) generators 
between the states in 3. The non-zero matrix elements 
are easy to find from the following equations giving the 
action of the generators on the states: 



Qs|o> = -||o>, 

■B|0> = ||0>, 



Os|l> = 0, 
BID = U) 



Q 3 |2> 

B\2) 



Q+|0) = |2), Q_|2> = |0>, 

V+\0) = \1), V-\2) = -\l) 



W+\l) = |2), 



W_|l> = |0>. 



(A6) 



These equations give us the matrices of the generators of 
sl(2|l) in the fundamental representation (for a generator 
G the matrix elements GV,-, i,j = 1,2,3, are defined by 



G\i) = j: j G ji \j)y. 



B = 



Q+ = 



v + = 



V- 




n = 



Q- 



W- = 



W+ = 




(A7) 



In Fig. g we give the weight diagrams for the funda- 
mental and adjoint representations which are useful for 
understanding some of our arguments. 

For the down-links the construction is similar. The 
oscillator realization of the sl(2|l) generators is now 



Q3 


= \if1fl + 


fVi + 1 )- 


Q+ 


= fJv 


q- = fin 


B 


= 




v+ 






V- 







w+ 



V2 



V2 



(b\fl-b}fl)- 



t Fh 



(A8) 



These operators satisfy the same commutation rela- 
tions as the ones on the up-links and act in the three- 
dimensional space spanned by the SU(2) singlets 



11} 

1 2) 



I vacuum; 
-W-\0) ■■ 

-Q-\o) -- 



V2 



(b\fl-bl^)\o), 



(A9) 



These singlets form the representation 3 of the sl(2|l) 
algebra dual to the fundamental 3. Note that the state 
|1) contains odd number of fermions, and, therefore, has 
negative square norm: 

(1|1) = -1. (A10) 
The action of the generators on the states in the rep- 
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B 




FIG. 2: Weight diagrams of sl(2|l). We show two doublets and the adjoint of the subalgebra gl(l|l) in the adjoint representation 
diagram. 



rcsentation 3 is easily found to be 

4io) 



Q 3 |6> 

B\Q) 
0+]2) 

W+\2) 



2 

-|0), 

|5>, 
-11), 



Os|l) =0, 

s|i) = -|I), 

Q-lo> = -|2>, 



Q 3 |2)=--|2), 



B\2) 



V-\l) 
W_|0) 



■|2), 



-5l2>, 



(All) 



This gives the matrices for the generators in 3: 



B = 




Q 3 = 



<2- 



w_ = 



w+ 




(A12) 



APPENDIX B: BOUNDARY SUSY 



We demonstrate in this appendix that the introduction 
of a reflecting boundary preserves the full sl(2|l) SUSY. 



First we note some useful relations satisfied by the 
bosons and fermions defined in section [II C| . For any 
function F, all bosons and fermions (denoted by c, c') 
except the negative norm ones satisfy the commutation 
relations, 



[c,:F(c\c):]=:^F(J,c):, 



[c\:F{c\c): 



:F{c\c) 



(Bl) 



where the : : denotes normal ordering. The negative 
norm operators satisfy, 







[c,:F(c\c):] = -:—F(c\c):, 
ac< 



[d,:F(J,c)-] =:F(c f ,c) 



dc 



(B2) 



One can first w^ite the transfer matrix for a single A 
node in the bulkED. Since the scattering at the node is 
diagonal in spin indices (see Eq. fll^)), the node transfer 
matrices are products of two independent transfer matri- 
ces for each spin direction: 



n t *° 

<T=U 



Ta]T A i. 



(B3) 



where 



exp (t Aa {ftfl + blll)){l-t A S* 

X exp (-t A a(f<rfa + b a b a )) , 

n f(J + n ba + nj a + n hr7 . 



(B4) 
(B5) 



16 



Let us also introduce the following notation: 

T+ = H exp {t Aa {flft + ftt gt )) , (B6) 

T = n(l-*L) in % (B7) 

T_ = JJcxp (-tArifaf* + b*ba)) , (B8) 

so that T A = T+T T_. 

The three terms correspond respectively to the cre- 
ation, propagation and destruction of boson and fcrmions 
on evolution along the vertical direction. Similar ex- 
pressions can also be written for the B nodes. In the 
spin-rotation invariant case for any particular realization 
of the disorder i n th e sca ttering matrices, using the re- 
lations in Eqs. (Bl, B2), it can be checked that each 
node transfer matrix in Eq. (B3) commutes with the 
sum of the eight generators of the superalgebra sl(2|l) = 
osp(2|2) (see Appendix defined on the up link and 
down link on which the node transfer matrix acts. 

Having defined the bulk nodes, we now consider the 
network on a semi-infinite half-plane with a fully re- 
flecting boundary cither along the horizontal direction 
or the vertical direction. Although it is clear that phys- 
ical quantities cannot depend on whether the boundary 
is defined along the horizontal or the vertical direction, 
the two cases have to be studied very differently within 
the second-quantized formalism. This is because of the 
fact that we have singled out the vertical direction as 
time and the tensor product of Fock spaces on which U A 
and Ub act is defined along a particular horizontal row of 
links. For definitcness, let us assume that the boundary 
is always composed of A nodes. 

We first consider the case of a reflecting boundary 
along the vertical (time) direction as shown in Fig. [y. 
In this case, we retain the periodic boundary condition 
along the time direction and hence also the supertracc 
STr in the correlation functions. Only the node transfer 
matrices on the boundary have to be changed to account 
for the complete reflection at the boundary. This can be 



implemented by setting t A = in Eq. (B4). In this case, 
the boundary node transfer matrix reduces to the triv- 
ial identity operator and consequently the operators U A 
and Ub still commute with all generators of the sl(2|l) 
superalgebra. 

As mentioned before, we could have equivalcntly cho- 
sen the boundary to be along the horizontal space direc- 
tion, that is along a single time slice. In this case, we 
will have to first replace the supertrace STr by the ma- 
trix element w.r.t the global vacuum state |0). Next we 
will have to modify all the node transfer matrices along 
the boundary by setting t A = 1 in Eq. (B4) and also 
consider only T + since no bosons or fcrmions can be cre- 
ated or propagated across the boundary. Note that the 
90' rotation of the boundary changes the corresponding 
tA- Hence the single spin single node transfer matrix at 
the boundary is: 



T A = Y[cxp[(p a p a + blbi)} 



(B9) 



This operator commutes only with the four elements, 
B + B, Q 3 + Q 3 , W+ + W+ and V- + VL which form 
the subalgebra gl(l|l). This seems to contradict the pre- 
vious observation that a reflecting boundary along the 
vertical direction preserves the full sl(2|l) SUSY. This is 
reconciled by the fact that in the former case, we took 
the supertrace STr with a trivial boundary node transfer 
matrix while here we need to consider the action of the 
boundary transfer matrix T A in Eq. ( p39|) on th e global 
vacuum |0). Using the relations in Eqs. ([Bl, B2), we can 
check that the state ?a|0) is a singlet under the action 
of the sl(2|l) symmetry on the two links involved, that 
is, it is annihilated by the sum over two links of all eight 
generators of the sl(2|l) superalgebra. Thus the full su- 
persymmetry is restored within a lattice spacing from the 
boundary and the result matches with the previous case. 
For simplicity, in the main text, we always assume that 
the boundary is along the vertical direction as shown in 
Fig. [|. This enables us to retain the global supertrace 
STr in all the expressions. 
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